home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / dlaed5.z / dlaed5
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. DDDDLLLLAAAAEEEEDDDD5555((((3333FFFF))))                                                          DDDDLLLLAAAAEEEEDDDD5555((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DLAED5 - subroutine computes the I-th eigenvalue of a symmetric rank-one
  10.      modification of a 2-by-2 diagonal matrix   diag( D ) + RHO  The diagonal
  11.      elements in the array D are assumed to satisfy   D(i) < D(j) for i < j
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE DLAED5( I, D, Z, DELTA, RHO, DLAM )
  15.  
  16.          INTEGER        I
  17.  
  18.          DOUBLE         PRECISION DLAM, RHO
  19.  
  20.          DOUBLE         PRECISION D( 2 ), DELTA( 2 ), Z( 2 )
  21.  
  22. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  23.      This subroutine computes the I-th eigenvalue of a symmetric rank-one
  24.      modification of a 2-by-2 diagonal matrix
  25.  
  26.      We also assume RHO > 0 and that the Euclidean norm of the vector Z is
  27.      one.
  28.  
  29.  
  30. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  31.      I      (input) INTEGER
  32.             The index of the eigenvalue to be computed.  I = 1 or I = 2.
  33.  
  34.      D      (input) DOUBLE PRECISION array, dimension (2)
  35.             The original eigenvalues.  We assume D(1) < D(2).
  36.  
  37.      Z      (input) DOUBLE PRECISION array, dimension (2)
  38.             The components of the updating vector.
  39.  
  40.      DELTA  (output) DOUBLE PRECISION array, dimension (2)
  41.             The vector DELTA contains the information necessary to construct
  42.             the eigenvectors.
  43.  
  44.      RHO    (input) DOUBLE PRECISION
  45.             The scalar in the symmetric updating formula.
  46.  
  47.      DLAM   (output) DOUBLE PRECISION
  48.             The computed lambda_I, the I-th updated eigenvalue.
  49.  
  50.  
  51.  
  52.  
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.